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The theoretical demonstration of spin transfer torque in metallic spin valves (SV) ten 
years ago [1] gave a new breath to giant magnetoresistance related studies [2], promising 
exciting new applications in non-volatile memories technology [3] and radio-frequency 
oscillators [4]. A number of fundamental studies in metallic spin valves revealed the 
different properties of spin torque and led to a deep understanding of current-induced 
magnetization dynamics [5, 6, 7, 8, 9]. Particularly, several theoretical studies described 
the structure of the torque in metallic magnetic multilayers and showed the important 
role of averaging due to quantum interference, spin diffusion and spin accumulation 
[10, 11, 12]. 

Since the first experimental evidence of spin-dependent tunnelling [13], magnetic 
tunnel junctions (MTJs) have attracted much attention because of the possibility to 
obtain large tunnelling magnetoresistance (TMR) at room temperature [14]. The 
possibility to use MTJs as sensing elements in magnetoresistive heads, as non-volatile 
memory elements or in reprogrammable logic gates has also stimulated a lot of 
technological developments aiming at the optimization of MTJs' transport properties 
and their implementation in silicon-based circuitry [3, 15]. Because of these applications, 
MTJs have been intensively studied and the role of interfaces [16], barrier [17], disorder 
[18] and impurities [19] have been addressed in many publications [21]. The recent 
achievement of current-induced magnetic excitations and reversal in MTJs [20] has 
renewed the already very important interest of the scientific community in MTJs. 

The observation of spin transfer torque in low RA (resistance area product) MTJ 
using amorphous [20] or cristalline barriers [15, 22] opened new questions about the 
transport mechanism in MTJ with non collinear magnetizations orientations. As a 
matter of fact, whereas the current-perpendicular-to-plane (CPP) transport in SV is 
mostly diffusive and governed by spin accumulation and relaxation phenomena [11, 12], 
spin transport in magnetic tunnel junctions is mainly ballistic and governed by the 
coupling between spin-dependent interfacial densities of states: all the potential drop 
occurs within the tunnel barrier. J. C. Slonczewski first proposed a free electron 
model of spin transport in a MTJ with an amorphous barrier [24], deriving TMR, spin 
transfer torque (STT) and zero bias interlayer exchange coupling (lEC). This first model 
only considered electrons at Fermi energy, neglecting all non-linear tunnel behaviour 
(consequently, current-induced lEC was found to be zero). More recently, the author 
proposed a more general model based on Bardeen's Transfer Matrix (BTM) method [25]. 
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Another group presented at tight-binding model (TB) of a MTJ, giving more realistic 
band structures than the usual free electron model [26, 27]. These studies showed that 
spin torque should present an important bias asymmetry and the dissipative part of lEC 
(also called current-induced effective field) should be of the same order of magnitude 
than STT with a quadratic dependence on the bias voltage [26]. Finally, we note that 
in the same spirit as Ref. [28], Levy and Fert studied the role of hot electrons-induced 
magnons on STT in MTJ [29]. In recent experiments, the important relative amplitude 
of current-induced effective field compared to the spin torque term has been verified 
[30, 31, 32] but the role of magnons is still under investigation (in the first experiment 
the current-induced magnetization reversal occured while the TMR was quenched by 
magnons emissions [20]). These specific features show that tunnelling transport has a 
strong influence on spin transfer torque characteristics. 

We recently presented [12] a description of spin-dependent transport in a MTJ 
treated in a free-electron assumption, based on Keldysh non-equilibrium technique [33] 
applied to a MTJs with an amorphous barrier (such as AlOj,). This method is close to 
Ref. [24], although more general since we consider the contribution of all electrons lying 
below the Fermi energy. In the present article, we focus on the anatomy of spin transfer 
torque in such a MTJ, paying attention to the origin of the speciflc characteristics of 
this torque in the particular case of MTJ. The paper is organized as follows. In section 
2, we give a reminder of the origin of spin transfer torque, and the way to calculate it. 
In section 3, the formulation of spin-dependent currents and torques in non-equilibrium 
Green function formalism (Keldysh formalism) is developed. Section 4 presents the 
results of the model and describe the anatomy of spin torque in a magnetic tunnel 
junction, underlying the role of tunnelling process. 

2. Current-induced torques 

All along this paper, we consider the s-d model in which 5-electrons are itinerant and 
(/-electrons are localized and give rise to the local magnetization of the ferromagnet. We 
also assume that the d local moments remain stationnary. This model applies to the 
electron structures of ferromagnetic electrodes whose compositions lie on the negative 
slope side of the Slater-Neel-Pauling curve [34] (Ni, Co, NiFe, CoFe). No spin flip is taken 
into account. In a magnetic tunnel junction composed of two semi-inflnite ferromagnets 
separated by a tunnel barrier (see Fig. 1), majority spins and minority spins refer to the 
electron spin projection in the left ferromagnet respectively parallel or antiparallel to the 
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local magnetization. In this framework, the motion of 5-Iike electrons in a ferromagnet 
is represented by the non-relativistic Hamiltonian including s-d coupling: 

2 

H=^ + f/(r) + Jsd (^.^d) (1) 
2m V / 

where the first and second terms are the kinetic and potential energies, the third term 

is the s-d exchange energy, S^^ being a unit vector collinear to the local magnetization 

due to the localized d-electrons, Jgd the s-d exchange constant and is the vector of 

Pauli matrices in spin space. After some algebra [12, 35], it is possible to derive the 

equation of continuity of the spin density : 

where \I/ = (\E'l^, \E'-'^) is an arbitrary 2-dimension Hartree-Fock wavefunction. The two 
dimensions refer to majority (t)and minority (J.) spin projections of the Hartree-Fock 
wavefunction. Here, ~^ (r, t) refers to the local spin density (namely the local out-of- 
equilibrium magnetization due to the itinerant polarized electrons): 

r{r,t) = ^*{r,t)^-^^{r,t) (3) 

Defining 

= -— a'® V.*(r,t)} (4) 

where J'^ refers to the spin-current density, we obtain, in steady states: 

2 T 

V,J^(r,t) = -^5,x^(r,t) (5) 

Eq. 5 implies that the spatial transfer of spin momentum from the itinerant s- 
electrons to the localized c?-electrons (left-hand side of Eq. 5) is equivalent to a torque 
exerted by the transverse spin density on the local magnetization (right-hand side of 
Eq. 5). This equivalence has been demonstrated by Kalitsov et al. [27] in magnetic 
tunnel junctions using Keldysh formalism and TB description. 

In the following, we calculate spin transfer torque from the torque exerted by the 
transverse spin density on the local magnetization. The spirit of our calculation is 
depicted in the top panel of Fig. 1. The out-of-equilibrium magnetic tunnel junction 
is modelled by a "conductor" (in the sense that the tunnel barrier is not infinite) 
linking two magnetic reservoirs {Fl and Fr) with non collinear magnetizations and with 
different chemical potentials fiL and fiR [36] (/iL > /ij?)- A bias voltage V = {jiL — /U/?)/e 
is applied across this "conductor". One should consider all electrons with majority 
spins (solid arrows) and minority spins (dotted arrows), originated from left (rightward 
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Figure 1. Schematics of the magnetic tunnel junction with non colUnear 
magnetizations orientation. Top panel: spin-dependent out-of-equilibrium transport 
in a conductor linking two reservoirs Fl and Fr (whose electrochemical potentials 
are respectively and ^b.) with non collinear magnetizations orientations. The solid 
arrows represent the majority spins and the dotted arrows represent the minority 
spins. Middle panel: MTJ with non collinear magnetization orientations. Bottom 
panel: Corresponding energy profile of the MTJ. In free-electron approximation, the 
local density of states are parabolic for majority (solid line) and minority (dotted 
line) electrons with a splitting between the two spin subbands equals to the exchange 
interaction Jsd- 

arrows) and right electrodes (leftward arrows). In low bias limit (/il ~ /i/j), the charge 
transport can be approximately determined by the electrons originated only from the 
left electrode at the Fermi energy. 

In our case (middle panel of Fig. 1), the magnetic tunnel junction is composed 
of two ferromagnetic layers, Fl and Fr (made of the same material, for simplicity), 
respectively connected to the left and right reservoirs and separated by an amorphous 
tunnel barrier. The a^-axis is perpendicular to the plane of the layers and the 
magnetization of Fl is oriented following z: Ml = Ml^- The magnetization A?R of Fr 
is in the {x,z) plane and tilted from il?^ by an angle 6. In this configuration, the spin 
density in a ferromagnetic layer possesses three components : rft = {rrix, ruy, rriz). In Fl 
(we obtain the same results considering Fr), the transverse components are m^. =< > 
and niy =< >, where cr* are the Pauli spin matrices and <> denotes averaging 
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over orbital states and spin states, i.e. averaging over electrons energy transverse 
momentum it and spin states. Thus, the transverse spin density exerts a torque on 
the background magnetization a}l following two axes: 



R 



(6) 



One should introduce the previous formula in the usual Landau-Lifshitz-Gilbert equation 
to describe the modified dynamics of the magnetization A^l- 



dt dt 

where a is the Gilbert damping, 7 is the gyromagnetic ratio and /is is the Bohr 
magnetron. The two terms in the right hand side of Eq. 6 stand for two types of 
torques: —rriyKi l x l x A?^^ is the usual STT term (also called in-plane or parallel 
torque[26]) whereas m^Mi^ Jiji is the current-driven interlayer exchange coupling (also 
called field-like torque, out-of-plane or perpendicular torque[26]). The former vanishes 
at zero bias, whereas the latter exists even without current [24, 26, 27]. An explanation 
of the physical nature and origin of these two terms will be given in section 4. The 
transverse spin density in the left layer is then given by < 0"+ >=< cr"^ + ia^ > : 

+ iniy =< cr+ >=< ^ j ^ ^ ^ ^ ^ >= 2 < > (8) 

In other words, STT is given by the imaginary part of < cr"^ >, while lEC is given 
by its real part. One can understand the product < \Ef*1^\Ef-'^ > as a correlation 
function between the two projections of the spin of the impinging electron. In ballistic 
regime, an electron impinging on a ferromagnet with a spin polarization tilted from the 
background magnetization will precess around this magnetization [10, 27]. Locally, its 
two projections j and | following the quantization axis (defined by the background 
magnetization) will be non-zero. Then, the electron will contribute locally to the 
transverse spin density rrix and rriy. If the electron spin is fully polarized parallel or 
antiparallel to this magnetization, no precession will occur and its contribution to the 
transverse spin density will be zero. 

We remind that we defined majority (minority) states as the spin projection parallel 
(antiparallel) to the magnetization of the left electrode. Then, < \E'*Tvl>i > will be the 
fraction of electrons whose spin is following x (real part) and y (imaginary part) in spin 
space. 
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3.1. Keldysh Green functions 

As explained previously, in Keldysh out-of-equilibrium formalism [33, 36], any physical 
quantity should be calculated considering the contribution of the electrons originated 
from the left reservoir and from the right reservoir (top panel of Fig. 1). Then, an 
out-of-equilibrium Green function G{r, t, r', t') (or Keldysh Green function) is defined 
as a superposition of these two contributions: 

G (r, t, r', t') = fL^L (r, t) (r', t') + fn^R (r, t) (r', t') (9) 

where ^'L(ij) (^7 1) are the electron wavefunctions originated from the left (right) reservoir 
at the location r and time t and are the Fermi distribution functions in the left 

and right reservoirs. 

Thus, the Schrodinger equation of the magnetic tunnel junction is: 

where the vector in Pauli matrices space : 0^ = (cr^', a^, cr^)^, E is the electron energy, 
U is the spin-independent potential along the junction: 

Jsd (^'o^-Sd^ = JsdO'^ and U = Ep for x < xi 

^\ X — X\ 

~^ .Sd] = and Uix) = Uq eV for Xi < x < X2 

J X2 — Xi 

Jsd (^-Sdj = Jsd (<7^ cos6' + sin 9) and U = Ep — eV for x > X2 



We consider that the potential drop occurs essentially within the barrier and we 
apply a low bias voltage compared to the barrier height (V « U/e). This allows to use 
WKB approximation to determine the wavefunctions inside the barrier. Furthermore, 
the free electron approximation implies parabolic dispersion laws which also restricts 
our study to low bias voltage. 

To describe the spin-dependent transport within the MTJ, we define the 
wavefunctions \E'°^ *^'^''(r,e), where e = Ep — E and E is the tunnelling electron energy. 

'•'^■'(r, e)p is the probability that an electron originated from electrode i, at the 
energy e, initially in spin state a, possesses a spin projection a' at the location r. For 
example, an electron initially in majority state, originated from Fl, is described by six 
wavefunctions along the structure: 
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Vh 

in the left electrode Fl {x < xi), 

^i(T) ^ c'^^(a:i,a:) + 4-^(a:,Xi) 
^ g(x) 



where E{xi, Xj) = exp ^ J^' q{x)dx^ , in the tunnel barrier (xi < x < X2), 



in the right electrode Fr (x > X2). A;i(2) and A;3(4) are the wavevectors for majority 
(minority) spin projection in the left and right electrodes, whereas q{x) is the spin- 
independent wavevector inside the tunnel barrier. Connecting the wavefunctions and 
their derivatives at the interfaces, we obtain the 24 wavefunctions (two spin projections 
and two reservoirs). These wavefunctions are given in Appendix. 

In the 2-dimensionnal Hartree-Fock representation, spin-dependent current and spin 
density are defined using the out-of-equilibrium lesser Keldysh Green function: 



G-J{r,r')= de I 



a'(T)* 



R 



(11) 



where //, = /°(e) and /r = /°(e + eV), and /°(e) is the Fermi distribution at 
K. For conveniency, we use the mixed-coordinate system (x, 7?), where it is the 
momentum parallel to the plane and x is the coordinate perpendicular to the plane. 
With r = (x, 7?), we get: 



G-J (r, r') 



ao 



G^^,{x,x')dlt 



(12) 



where ao is the lattice parameter of the electrodes [37]. Spin transfer torque (STT, T\\) 
and interlayer exchange coupling (lEC, T^) can now be determined from Eq. 8, whereas 
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+ = — < a+ >= 2—-^ [ [ 6-7+ (x, x, e)KdKde (13) 
AtB fJ'B (27r)^ J J '-^ 

^ j j [G'^^ix.x.e] - G~^^{x,x,e)\KdKde (14) 



he 




d___d_ 

dx dx' 



G^^~^^^^-^{x,x',e)\j:=x'KdKde (15) 



J=J^ + Jl (16) 

G|^^(x, x,e) and Gj^{x,x,e) are the energy-resolved local density-of-states (LDOS) for 
up- and down-spins respectively, whereas / G^.^{x,x,e)de and / Gj^{x,x,e)de give the 
number of up- and down-electrons at the location x along the structure. 

3.2. Calculation of spin transfer torque 

As demonstrated in Eq. 5, it is possible to calculate spin transfer torque from the 
divergency of spin current density or from the spin density itself. We now demonstrate 
that this relation holds in our model. Spin current densities and spin density are defined 
as [10]: 



my = —I 



^.^ = -^^{**'^-**'^} (20) 



— * 

Im ox ox 

We evaluate these quantities for electrons originating from the left reservoir in the 
left electrode (x < xi). The equations are given in Appendix. The spin densities for 
majority (j) and minority (|) electrons are: 

(^-i{kx+k2){x-xx) _ ^*'\ ^i{kx-k2){x-xx) \ 
d^ ^^'^J ^^^^ 

(^-i(k\+k2){x-xx) _ ^-i{ki-k2){x-xx) \ 
d^ ^ ] ^^^^ 

len I ^'^'^^ 

I Q-i{ki+k2){x-xi) _ y,*l^~i{ki~k2){x-xi) \ 

= —8iqiq2{k3 — k^) sin 6 ( j-^ c.c. ) (24) 
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m], + = 8gig2(^3 — ^4) sin 9 



(25) 



X 2 



-i{ki+k2)(x-xi) 



den 

8zgig2(^3 — ki) sin 6* 



c.c. 



r 



1 , '1 



+ 



T 



den den* 



den* den 



-i(ki~k2){x—xi) 



C.C. 



(26) 



By the same way, we evaluate the spin current density for majority and minority 
spins: 

J^2 I ^-i(ki+k2){x-xi) ^i{ki-k2){x-xi) 

Jf" = -8qiq2-z—(k3 - sin6' -ik2 ; h ik2— ; 

2m \ den den 



^i{ki+k2)(x~xi) ^Tg-i(fci-fc2){a;-xi) 

+iki ; h iki— — 



den* 



den* 



(27) 



^2 / gi(fci+fc2)(x-a;i) j.lQi{ki-k2){x-xi) 

Jx^ = -8qiq2 — ik3 - ki) sin^ ik2 ^— ^ + ik2- ^— ^ 

2m \ den den 



g-j(fci+fc2)(x-a;i) ^-i{k\-k2){x-xi) 

—iki ; h iki — 



den 



den 



(28) 



J^^ = -Sqiq2^{k3 - ki) sin 6* I ik2- 



'2m 



+iki- 



-i{ki+k2){x-xi) 

den 

,i(ki+k2){x—xi) 



-lk2^ 



T i(A:i-fc2)(x-xi) 



r, e 



den* 



+ ^/cl — 



den 

^-i{k\-k2){x-x\) 



den* 



(29) 



^2 / gi(fci+A:2)(2:-xi) ^-Lgi(fci-fc2)(a;-xi) 

J'-^ = -8gig2 — (^3 - ^4) sin^ -z/cg 1— ^ «A;2^ 1— ^ 

^ 2m \ den* den 

^-i{ki+k2){x-x-i) ^-i{ki-k2){x-xi) 

—iki ; h iki— — 



den 



den 



(30) 



Taking the imaginary (and real) part of the right-hand-side of the above equations, 
we obtain, similarly to Eqs. 25 and 26: 



1 Aki + k2) . . 
Jx = -8qiq2^{k3 - k^) smO 



+ 



den den* 



^i{ki-k2)(x-xi) _|_ ^^^^ 



(31) 
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2m 



^i{k-i+k2){x-x\) 

den* 



CO. 



{k, - k2) 
4 



den den* 



gj(fci-fc2)(a:-xi) _ ^^^^ 



(32) 



The divergency then gives: 



dx 



fe2 l2 _ 1,2 

-8iqiq2 — {k3-k^) sm9 

Zm 2 



den den* 



J{ki-k2){x-xi) _ 



C.C 



dx 



^2 _ ^2 
1^1 1^2 



8qiq2 — {k^-ki)sm.6 

2m 2 



-i{ki+k2){x-xi) 



den* 



+ C.C. 



+ 



den den* 



^i{ki-k2){x-xi) _|_ 



(34) 



Setting Jsd 



2m 2 ' 



Eqs. 25, 26, 33 and 34 give the following relation: 



Ox ^ sd^y 

ox 

=^ VJ' = JsdM X m (35) 

Then, the relation 5 can be derived analytically in the free electron approach. This 
relation does not depend on the particular description adopted (Tight-binding or free 
electron approximation) but emerges from the definition of the considered Hamiltonian 
itself. 



4. Results and discussion 

To illustrate the above calculation, we use material parameters adapted to the case 
of C0/AI2O3/C0 structure: the Fermi wavevectors for majority and minority spins 
are respectively kp = 1.1 A~^, kp = 0.6 A~^, the barrier height is U — Ep = 1.6 
eV, the effective electron mass within the insulator is meff=0.4 [38] and the barrier 
thickness is d=0.6 nm. These parameters have been choosen to fit the experimental I-V 
characteristics of the magnetic tunnel junctions studied in Ref. [31]. In all this section, 
the magnetizations form an angle 6=90° between them. We will justify this choice in 
the following. 
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4-1- Anatomy of spin transport 

Although spin-dependent tunnelling is a well known process, the description we give 
here is of great importance to understand the specific characteristics of spin transfer 
torques in tunnelling transport. In this part, we will consider the linear approximation 
in which the bias voltage H is low enough so that the current is due to Fermi electrons 
injected from the left electrode. When the electrodes magnetizations are non coUinear, 
the electrons are no more described as pure spin states, but as a mixing between majority 
and minority states. For example, let us consider one electron from the left reservoir, 
initially in majority spin state, impinging on the right electrode (see Fig. 2 - step 1). 
The first refiection (step 2) at the Fl/I interface do not introduce any mixing since 
the insulator is non magnetic. However, when (the transmitted part of) this electron is 
reflected or transmitted by the second interface I/Fr (step 3), the resulting state in the 
right electrode is a mixing between majority and minority states since the quantization 
axis in the right electrode is different from the quantization axis in the left electrode. 
Then, the transmitted spin is reoriented and precesses (step 4) around the magnetization 
of the right electrode. Furthermore, the reflected electron (step 5) is also in a mixed 
spin state and precesses around the left electrode magnetization. In other words, after 
transport through the barrier, the electron spin is reflected/transmitted with an angle. 
This reorientation gives rise to spin transfer torque. 

Note that there is not reason why the electron spin should remain in the plane of 
the electrodes magnetizations. We will see that after the reorientation, the electron spin 
possesses three components in spin space (and so two transverse components). 

4.1.1. Tunnelling transport We are first interested in the spin-dependent reflectivity 
^CT(cr') g^j^^ transmittivity T"'^'^'^ for electrons at the Fermi energy from the left electrode. 
Let us consider an electron initially in majority spin state (|). Its wavefunction will be 
described by a plane wave in the left electrode : 



The mixing between majority and minority spin states can be expressed through 
mixing reflectivities R)^ and R'^^ and transmitivities T^^ et T^^ , so that: 



,ik\{x—x\) 




where: 
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Figure 2. Schematics of the principle of spin transport in a magnetic trilayer with non 
colUnear electrodes magnetizations. Step 1: the electron spin is polarized along the 
magnetization of the left electrode. Step 2: After the first reflection/transmission by 
Fl/I interface the reflected and transmitted parts remain in a pure spin state. Step 
3: The reflection/transmission by the second interface I/Fn reorientes the electron 
spin. Step 4 and 5: The transmitted and reflected spins precess around the local 
magnetization. 



R^^ = \rl\^ (36) 

= lQ\^_lSl^^l^,i^0\2 (37) 

m^j^jden 

AT/*T(T) ^,T,T(T) 
Tn = |vl;m)^!:z__g,*T(T)^| (38) 

^ ^ dx ^ dx ^ ^ 

^ ^ dx ^ dx ^ ^ ' 

^ is evaluated in the right electrode and given in Appendix. By the same 
way, we can define the transmittivity and reflectivity of an electron initially in minority 
spin state. Fig. 3 displays the K-dependence of R and T (we omit the superscripts for 
simplicity), where k is the wavevector component in the plane of the layers. 

More than 97% of the majority and minority spins are reflected conserving 
their spin projection, whereas less than 3% are transmitted without spin flip. This 
reflectivity (transmittivity) reaches a minimum (maximum) at perpendicular incidence 
and increases (decreases) quickly with k. Note that T^^ and T^^ are equal due to 
the particular configuration of the electrodes magnetizations (^=90°). Thus, after 
interaction with the barrier, only a very small part of the spin is flipped (the flipped 
spins have to tunnel through the barrier twice) : less than 2.7 x 10^'^% of the reflected 
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K(a-i) 

Figure 3. Reflectivity (top panel) and transmittivity (bottom panel) as a function of k. 
The solid and dotted lines represent the spin conserving reflectivity and transmittivity 
for initially majority and minority spins respectively (left axes); the dashed and dotted- 
dashed lines represent the mixing reflectivity and transmittivity for initially majority 
and minority spins respectively (right axes). The applied bias voltage is Vf, = 0.1 V 
and 61=90°. 

wave has flipped its initial spin. 1.6x10"^% of the electron spins initially in minority 
states reverses its spin during reflection. 

Thus only a very small part of the injected polarized wave is flipped during the 
tunneling process. However, this does not mean that spin transfer torque is small in 
MTJ; as a matter of fact, only coherent mixing states will contribute to transverse spin 
density, generating spin transfer torque. 

Finally, we note that only electrons close to the perpendicular incidence contribute 
significantly to the current. This has important consequences on the impact of quantum 
interferences on spin transfer. 
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4.1.2. Spin density and spin transfer torque In the linear regime under consideration, 
the three components of spin density in the left electrode can be described as follows: 



= Vl/Ta)vl;:i(i) 



_ 8qiq2{h - h) sin 6 

8qiq2{h - h) sinO 
mljjrden 



^i{ki+k2)[x-xi) _ ^'\^^-i{ki-k2)ix-xi) 



g-j(fci+/c2)(x-a;i) _ ^l* ^~i(ki~k2)(x~xi) 



(40) 
(41) 



m 



zL 



^t(T)^*T(T)_^|(T)^;i(T) 



(42) 



,Ta),Tr*i(T) 



8qiq2\/ki{kz ~ ^4) sin 6' 



den 



1 



^2iki{x-xi) 



+ r{e 



^ -2iki(x-x-i_) 



^ia)^*i(i) 



(43) 



8qiq2\/h{h - h) sinO 



m 



eff 



den 



— f j-*^q'^^'^2{x-Xi) _|_ ^j,^ 



-2ik2{x—xi) 



Observing m'_;l' in Eq. 40, one can distinguish two components: the first one 
is proportional to e±«('=i+'^2)(a;-xi)^ ^^^^ interference between the incident 

wave with majority (resp. minority) spin and the reflected wave with minority (resp. 
majority) spin; the second one is proportional to e^'^^''^~'''^^^^~^^^ and due to the 
interference between the reflected waves with majority and minority spins. We note 
that the flrst components of m^^^ and m^j^ are complex conjugated so that their sum 
is real. Then, the interference between the incident wave with majority spin and the 
reflected wave with minority spin does not contribute to STT but only to lEC. STT is 
then generated by the coherent interferences between reflected electrons with opposite 
spin projection (oc Q~^i>'i-f'2)ix-xi)y 

Concerning nizLj it is composed of one component proportionnal to e^*'^i(^~^i\ one 
component proportionnal to e^**^2(a:-a::i) g^j^^j constant as a function of x. The two 
formers are due to the interference between wavefunctions in the same spin projection 
but with opposite propagation direction while the latter is due to interference between 
wavefunctions in the same spin projection and the same propagation direction. 

Fig. 4 displays the details of the spin density components m^,, niy et (described 
in Eq. 40) in the left electrode as a function of x, when Vb = 0.1 V. possesses a quite 
complex behaviour with two periods of oscillation (the dashed lines show the enveloppe 
of the curve), whereas niy is reduced to a single oscillation (The oscillation period ki + k2 
vanishes when suming the contribution of majority and minority spins); oscillates 
around mean values represented by horizontal dashed lines. 
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Figure 4. Projections of spin density due to Fermi electrons in perpendicular incidence 
from the left electrode, as a function of the distance from the interface. Top panel: 
rrix component of spin density (solid line) ; the dashed lines are the enveloppes of the 
curve. Middle panel: my component of spin density. Bottom panel: to^ component 
of spin density due to initially majority (solid Hne) and minority (dotted Hne) spin 
projection; the dashed Hues are the mean values of the oscillations. The applied bias 
voltage is Vfc = 0.1 V. The vertical line on the right is the interface between the left 
electrode and the tunnel barrier. 

Note that the conservative part of lEC is only proprotionnal to e~^^'''^~^''^^^^~^'^\ But 
at non zero bias, the dissipative part of lEC is proportionnal to both e-«(fci+*^2)(x-a;i) ^^^^ 
^-t{ki-k2)(x-xi) ^ Then, at non zero bias, the electrons will not precess circularly around 
the background magnetization, but will present a more complex structure. 

Following the previous discussion about spin reorientation (see Fig. 2), it is possible 
to deduce the angles at which the electron spin is reflected by the barrier. We deflne 
the azimuthal angle azimuthal f] and the polar angle as indicated in the insert of Fig. 
5: 



Description of current- driven torques in magnetic tunnel junctions 



17 



7] = arctan — - (44) 

ml 

cp = arccos — (45) 

^J'mf:^ + vriy + m'f 

In definition of the vector ryf^, we only considered the coherent interferences between 
plane wave propagating in opposite direction (oc e^*'-'^^"^^-' and constant component 
of m^), as discussed above. Fig. 5 shows the ^-dependence of these angles at the 
interface Fi^jl (x=-3 A) for an electron spin initially in majority state and for different 
barrier thicknesses (top panel) and heights (bottom panel). The azimuthal angle r] varies 
between -64° to +77° while the polar angle remains very small (less than 0.2°), which 
means that the electron spin stays very close to the quantization discussed above. 

At K = 0.6 (corresponding to fc^), r/ = which indicates that the effective spin 
density lies in the plane of the magnetizations {^l-, ■ Finally, the polar angle does 
not vary with the distance, which means that the reflected electron spin precesses around 
Oz with a small angle 0. A "Bulk" spin transfer only occurs under the interferences of 
all the reflected electrons. 

The strong dependence of r/ as a function of the in-plane wavevector k, together 
with the dominant contribution of nearly perpendicularly incident electrons (see Fig. 3), 
implies that the effective electron spin, resulting from the averaging over all the incident 
electrons, possesses an important out-of-plane component. In other words, the effect of 
the spin-dependent tunneling is to strongly enhance the dissipative lEC component of 
the spin torque, compared to metallic spin valves. As a matter of fact, in SV the whole 
Fermi surface contributes to the spin transport so that the effective angle r] is very small 
[10]: the dissipative lEC is thus negligible. 

Note that increasing the thickness of the barrier only weak influence on t] and 
strongly decreases the amplitude of (the mixing reflection decreases since the barrier 
thickness increases, then reducing the transverse spin density). Furthermore, when 
increasing the barrier height, both amplitudes of the angles and r] decreases near the 
perpendicular incidence. These results are consistent with the reduction of spin mixing 
when decreasing the barrier transmittivity. 

Fig. 6 shows the dependence of the angles as a function of the s-d exchange 
constant Jgd for perpendicular incidence k = 0. Quite intuitively, the precession angle 
increases with Jsd whereas the initial azimuthal angle rj decreases in absolute value 
with Jsd- The spin-flltering effect (the selection between majority and minority spin 
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K(a-1) 

Figure 5. K-dependence of the reflection angles for an electron spin at Fermi energy, 
initially in majority spin state. Top panel: the barrier thickness is set to d=0.6 nm 
(solid line), d=0.8 nm (dashed line) and d=l nm (dotted line); U — Ep=1.6 eV. Bottom 
panel: the barrier height is set to U — Ep=1.6 eV (solid line), U — Ep=2 eV (dashed 
line) and U — Ep=3 eV (dotted line); d=0.6 nm. Insert: Deflnition of the angles (j) and 
77. The applied bias voltage is Vb = 0.1 V and 9=90°. 

during the reflection process) increases with Jgd so that ntc gets closer to the plane of 
the magnetizations. 

4.2. Spin Transfer Torques 

We now take into account all the electrons in the calculations (from the left and the 
right electrodes). Fig. 7 shows STT and lEC as a function of the angle 9 between 
the electrodes magnetizations, at H = and H = 0.1 V. It clearly appears that 
lEC and STT are proportionnal to sin^ (the deviation from sin 6* is minor than 10~^). 
This dependence is strongly different from what was predicted in metallic spin valves 
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Figure 6. Reflection angles as a function of the s-d exchange constant, for a Fermi 
electron initially in majority spin state. The parameters are the same as in Fig. 5. 

[11, 12, 39] and has been attributed [25] to the single-electron nature of tunneling. 
As a matter of fact, because of the important height of the tunnel barrier, all the 
potential drop occurs inside the insulator and spin accumulation (i.e. the feedback of 
the current-induced longitudinal spin density on the spin current) is negligible. In this 
case, the angular dependence of torque is determined by the angular dependence of the 
transmition matrix, as discussed in Ref. [25] and yields a sine shape. In the following, 
we will estimate the spin density for 6' = 7r/2. 

Note that, at zero bias, interlayer exchange coupling is still non-zero, contrary to 
spin transfer torque. The conservative part of lEC (lEC at zero bias) comes from the 
contribution of electrons located under the Fermi level. At zero bias, the currents from 
left and right electrodes are equal, but the electron propagation still corresponds to the 
scheme shown in Fig. 2: the mixing between majority and minority states induces a 
transverse component in the spin density. 

Fig. 8 displays the two components of transverse spin density as a function of 
the location in the left electrode. The interference process between polarized electrons 
yields a damped oscillation of lEC as presented in Fig. 8(a). We can distinguish two 
periods of oscillation Ti = 27r/ (^kp — kp^ and T2 = 27r/ (^k^ + kp^ whereas at zero bias, 
only T2 appears (see inset of Fig. 8(a)). This can be easily understood by considering 
electrons from left and right electrodes. Transverse spin density in the left electrode due 
to electrons from the right electrode is: 
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e (rad) 



Figure 7. Angular dependence of spin transfer (grey) and interlayer exchange coupling 
(black): conservative part (at zero bias - solid lines) and dissipative parts (bias 
dependent part - dashed lines). The dissipative parts are calculated at Vb — 0.1 V. 
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It is now possible to show that in the general expression of transverse spin density 



m 



+L 



m 



+L 



the terms proportional to e~'^^^^~^'^^'^^~^'^^ vanish at zero bias and m+ reduces to terms 
proportional to _ Furthermore, these last terms only give a real component 

since, as discussed above, the majority and minority components of rriy compensate 
each other. Consequently, at zero bias, only the conservative part of interlayer exchange 
coupling exists, due to the interference between incident and reflected electrons with 
opposite spin projection. But when the bias voltage is non zero, the transport becomes 
asymetric and the terms proportional to (Jq not compensate each other 

anymore and lead to two periods of oscillations as shown in Fig. 8(a). 

Spin transfer torque, proportional to my, only exits at non zero bias and possesses 
only one period of oscillation Ti (see Fig. 8(b)). It is worthy to note that the transverse 
components of spin density is damped by 50% within the first nanometers, and that 
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x(a) 

Figure 8. Total spin density as a function of the location in the left electrode: a) 
Current-induced interlayer exchange coupling - inset: Interlayer exchange coupling at 
zero bias voltage; b) Spin transfer torque. These quantities are calculated at Vf, = 0.1 
V. 

the amplitude of lEC is of the same order than STT. This decay lenght is very large 
compared to previous theoretical predictions [10, 39] and experimental investigations on 
SV [40]. As a matter of fact, the ballistic assumption holds for distance smaller than the 
mean free path (~ 5 nm in Co). In realistic devices, spin diffusion should increase the 
decay of STT and lEC. Another source of this difference compared to metallic SV is the 
fact that we consider perfect interfaces and no defaults in the barrier. First principle 
studies of realistic Co/Cu interfaces [41] showed that the mismatch of the electronic 
structure at the interface strongly reduces the transverse component of spin density. In 
MTJ, the non spherical nature of the spin-dependent Fermi surface [42, 4-3] should also 
dramatically alter the transverse spin density. This could also explain the fact that the 
amplitude of spin torque is two orders of magnitude higher than in experiments. 

Another characteristic specific to MTJ is that in our calculation we find that 
dissipative lEC is of the same order of magnitude than STT. This is coherent with 
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Figure 9. Out-of-equilibrium longitudinal spin density along the magnetic tunnel 
junction for majority (solid line) and minority (dotted line) electron spin projections. 
The bias voltage is 14 = 0.1 V. 



the theoretical results of Theodonnis et al. [26] as well as with the experimental studies 
of Petit et al. [31]. This can be attributed to the high K-selection due to the tunneling 
transport. We previously found that the contribution to torque strongly decrease with k 
(see Figs. 3 and 5) so that only electrons with small k strongly contribute to spin torque. 
In this case, the averaging of torques (and specifically lEC) will be less destructive 
than in metallic spin valves where all the Fermi surface is involved in the quantum 
interferences. 

Finally, Fig. 9 shows the out-of-equilibrium longitudinal spin density An defined 
as An^^^'> = n'^^^'>{Vb = 0.1) — ^^^'^■'^^(V;, = 0). An oscillates and asymptotically reaches a 
non zero value. This means that when the bias voltage is turned on, a non equilibrium 
spin accumulation builds up. However, this effective spin accumulation is very small 
{An^ — An^ ^ 10"^ electron/atom) and cannot infiuence spin current building. Then, 
neglecting the role of longitudinal spin accumulation (spin density) in MTJ is justified. 

4.3. Bias dependence 

The bias dependence of STT and lEC in MTJ also presents strong differences with SV. 
We first calculate the total spin torque exerted on the left electrode. Following the 
definition of Ref. [1] and Ref. [26], the total torque is: 

P — OO 

f total = / -wrdx = rixi) (50) 
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Figure 10. Bias dependence of interlayer exchange coupling (a) and spin transfer 
torque (b) for different values of s-d coupling: Jsd = 0.38 eV (open circles), Jsd = 0.76 
eV (filled circles), Jsd = 1.62 eV (open squares), Jsd = 2.29 eV (open triangles), 
Jsd = 2.97 eV (filled squares). Top inset: Bias dependence of STT for Jsd = 1-62 eV; 
the solid line was calculated following the usual way and the symbols were calculated 
using Eq. 51. 

Fig. 10 displays the total interlayer exchange coupling (a) and spin transfer torque 
(b) as a function of the applied bias voltage, for different values of the s-d exchange 
parameter J^^. Our results are consistent with Theodonnis et al. [26]. The dissipative 
lEC is quadratic whereas STT is a combination between linear and quadratic bias 
dependence. In Ref. [26], the authors proposed a general formula, derived from 
Slonczewski circuit theory [39], linking total spin transfer torque with interfacial spin 
current densities [26]: 

^1 = 2 ^ ' 

where Jap{p) interfacial spin current densities when the electrodes 

magnetizations are antiparallel and parallel respectively (see the definition in Ref. [26]). 
The authors claimed that this relation should hold for any electronic structure, so any 
transport description. As a matter of fact, the top inset of Fig. 10(b) shows STT 
calculated using Eq. 50 (solid line) and Eq. 51 (symbols). It shows very good agreement 
between the two members of Eq. 51. 

Experimental studies by Cornell's group [32, 44] demonstrated a linear variation of 
spin transfer torque as a function of the applied bias voltage. This linear variation is also 
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Figure 11. Bias dependence of interlayer exchange coupling (a) and spin transfer 
torque (b) for Jsd ~ 1-62 eV and different values integration depth: t = A(open 
squares), i = 4 A(filled triangles), t = 10 A(filled circles), t = oo A(open circles). 

usually assumed in interpreting excitations studies [30, 31]. Moreover, the very recent 
article of Sankey et al. [32] seems to confirm the fact that the dissipative exchange 
coupling is quadratic as a function of the bias voltage. Finally, note that a change of 
sign of spin transfer torque at high positive bias voltage is expected, consistently with 
Ref. [26]. The STT change of sign should be observed in MTJ with low enough barrier 
height and high breakdown voltage (MgO seems a good candidate). Nevertheless, more 
technological development is needed to fabricate such junctions. 

Eq. 50 assumes that all the transverse spin density has been absorbed in the free 
layer. However, in very thin free layer, one can expected that transverse spin density is 
not fully absorbed when leaving the free layer. In this case, one should consider that the 
free layer is finite. Fig. 11 displays the bias dependence of lEC and STT for different 
integration depths t (namely, different layer thicknesses): 

Xl — t 

-VJ'dx = r{xi)-r{xi-t) (52) 

1 

The dependence can change drastically and lEC can even change its sign (note that 
STT keeps its general shape). These dependences are strongly affected by the tunnel 
barrier characteristics and one should to be careful in the analysis of bias dependence. 



partial 



Description of current- driven torques in magnetic tunnel junctions 25 
4.4- From weak ferromagnetic to half-metallic tunnel junction 

To conclude this article, we study the dependence of the total spin transfer torque and 
interlayer exchange coupling as a function of the energy of the bottom of the minority 
electrons conduction band e^^, as indicated in Fig. 11. This energy is defined from the 
Fermi energy as: 

where is the absolute energy of the bottom of the conduction band. In the present 
study, we vary e-*^, keeping and Ep constant. When is close to e\ kp f» kp, the 
metallic electrodes loose their ferromagnetic nature. For ~ 0, the Fermi wavevector 
for minority electrons becomes smaller and the current polarization is strongly enhanced. 
In this case, we expect an important spin transfer torque. When > 0, kp becomes 
imaginary and the electrodes behave like a tunnel barrier for minority spins. Increasing 
e-*^ increases the evanescent decay of minority wavefunctions in the electrodes. Then, the 
product < \E'*1^\1/J^ > still exists so that spin torque is non zero and decrease exponentially 
from the interface. 

Fig. 12 shows the amplitude of total STT and current-induced lEC in the 
three different regimes: weak ferromagnetic electrodes (WFM), strong ferromagnetic 
electrodes (SFM) and half-metallic electrodes (HM). As expected, in ferromagnetic 
regime, STT and dissipative lEC increase until e^^ = (vertical line). When e^^ becomes 
positive, the bottom of the conduction band of minority electrons lies above the Fermi 
level: no minority electrons can propagate because only evanescent states exist near the 
interfaces for this spin projection. However, STT and dissipative lEC do not vanish but 
reach a plateau which slowly decreases to zero when increasing Jgd (not shown). 

To understand this behaviour, we calculated the spatial dependence of the 
transverse spin density in the free layer. Fig. 13 shows the transverse spin density in a 
usual ferromagnet, = —1.37 eV (which corresponds to Jsd = 1-62 eV), as a function of 
the distance from the interface in the left electrode. The oscillation possesses the same 
characterisics than discussed above and we observe that the transverse spin density is 
damped far from the interface. When decreasing e^, the interfacial spin density increases, 
due to strong spin filtering at the interface (strong spin-dependent selection), as shown 
on Fig. 14. 
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Figure 12. Spin transfer torque (solid line) and dissipative interlayer exchange 
coupling (dotted line) as a function of s-d exchange coupling. The vertical line shows 
the limit between ferromagnetic (weak ferromagnetic -WFM- and strong ferromagnetic 
-SFM-) regime and half-metallic regime. 




Figure 13. Transfer spin density (black line) as a function of the distance in the left 
ferromagnetic electrode in a usual ferromagnetic regime. We set = —1.37 eV and 
Vb = 0.1 V. 
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Figure 14. Transfer spin density (black line) as a function of the distance in the left 
ferromagnetic electrode in a strong ferromagnetic regime. We set = —0.38 eV and 
Vb = 0.1 V. 

But when e^^ changes its sign, only majority electrons can propagate and the 
transverse spin density is (see Eqs. 21-24): 

(g-i(fci+fc2)(x-x-i) _ „i{ki-k2){x-xi)\ 
^ ^^^^ 

(^-i{ki+k2){x-xi) _ *1 i{ki-k2){x-xi)\ 
) ^ ^^^^ 

Considering Fermi electrons at perpendicular indidence, very small bias voltage 
{eV ~ 0) and imaginary minority electron spin wavevector, /c2(4) = ik, we obtain 
straightforwardly: 



^-iki{x—xi) ^iki{x—xi) 



IGgigse'^^"-"^) sin 9 ^{{ks - ik) — ^ } (56) 



T 




-legigse^'^^-"^^) sin^ Q{{k3 - ik) — ^ } (57) 



The transverse spin density is a product between oscillating function of ki and 
exponentially decaying function of k. Fig. 15 shows the spatial evolution of the 
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Figure 15. Transfer spin density (black line) as a function of the distance in the left 
ferromagnetic electrode in half-metallic regime. We set = 19 eV and Vb = 0.1 V. 

transverse spin density in the case of a half-metallic tunnel junction. All the oscillations 
are damped very quickly so that the only important contribution to torque comes 
from the interface. Contrary to usual MTJ (where both bulk averaging due to spatial 
interferences and interfacial spin reorientation contribute to spin torque), in a strong 
half-metallic tunnel junction all the torque comes from spin reorientation due to spin- 
dependent reflection. In this last case, the contribution of the spatial averaging between 
all the impinging electrons (^-summation) is reduced compared to interfacial spin 
transfer. 

5. Conclusion 

A free-electron s-d model has been proposed to analyze spin transfer effects in magnetic 
tunnel junctions with amorphous barrier and non collinear electrode magnetizations. 
We first studied the anatomy of spin transport in such MTJ, showing that only a small 
part of the current undergoes spin-fiipping due to the non coUinearity of the electrode 
magnetizations. This corresponds to only a small deviation of the refiected spin from 
the local magnetization. Nevertheless, we showed that this small amount of precessing 
spin gives rise to an important transverse spin density leading to spin torque. 
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We also showed that the tunnel barrier acts like an incidence filter which increases 
the contribution of the electrons impinging with angle close to the perpendicular 
incidence. This K-selection is at the origin of an important lEC, contrary to what 
is observed experimentally in metallic spin valves. The ballistic transport dominating 
the tunnel transport in MT J is expected to induce large oscillations of STT and lEC as a 
function of the distance from the interface. If the oscillation period is large compared to 
the exchange length and one will observe a twist of the magnetization in the thickness of 
the layer. Otherwise, if the oscillation period is short compared to the exchange length, 
one will observe the torque integrated over the layer thickness. 

The bias dependence of spin transfer torque shows a strong asymmetry and a change 
of sign at positive bias voltage. This results is coherent with tight-binding calculations 
[26]. However, we saw that this model is strongly limited to small bias voltage because 
of the simplicity of the adopted band structure. 

Finally, we analyzed STT and lEC when varying the s-d exchange coupling and we 
demonstated that the torque still exists in MTJ composed of half-metallic electrodes, 
due to spin-dependent reflections. However, for inflnite half-metallic MTJ (for infinite 
s-d coupling), it is shown that STT and dissipative lEC vanishes to zero. 

Furthermore, several numerical studies have shown that, even in amorphous barrier, 
the interfaces composition and specially the presence of interfacial oxygen have a 
very deep influence on the spin polarization and thus on TMR and STT [42]. The 
recent development of MgO-based MTJ in spin transfer studies reduced the interest 
in amorphous barriers. However, amorphous barriers have the ability to present a 
simple physical framework which can constitute a basis to understand spin transfer in 
MTJ. Nevertheless, because of its more complex band structure and spin-filtering effect 
associated with the symmetry of wavefunctions, microscopic analysis of spin transfer 
in MgO-based MTJ would present exciting fundamental characteristics even on spin 
transfer effects [43]. 
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In this appendix, we give the analytical formulae for the spin dependent wavefunctions 
in the MTJ. Some functions which will be used in the description of this wavefunctions 
are first defined: 



E{xi,Xj) = exp / q{x)dx 

Jxi 

En = E{xi,X2) 

where Ep is the Fermi energy, U is the height of the barrier, V is the bias voltage and 
e = Ep — E, E being the energy of tunnelling electron. We define: 



^(gi, ki, q2, kj) = En{qi - iki){q2 - ikj) - (gi + iki){q2 + ikj) 

(f){qi, ki, q2, kj) = En{qi + iki){q2 - ikj) - E~^{qi - iki){q2 + ikj) 

den = ^(gi, ki,q2, /c3)*(gi, k2,q2, k4){l + cos6') + ^(gi, k2,q2, A;3)^(gi, h,q2, k^)^! - cos 6*) 
T 1 

ri = [(f){qi,ki,q2,k3)'^{qi,k2,q2,ki){l + cos6) + 0(gi, A;i, ga, /c4)^(gi, /c2, 92, A;3)(l - cos6')] 
^3 = — [0(g2, ks, gi, A;i)^(gi, /c2, 92, + cos 61) + 0(g2, k^, gi, A;2)^(gi, h, ^2, h){^ - cos9)] 

Electrons initially in the left electrode have the following wavefunctions along the 
structure : 



2m 



{U - Ep) 




gi = g(xi) 
92 = q{x2) 





CO < X < Xi) 




1 




(— oo < a; < xi) 



8gig2 {ks — /C4) sin 9 



den 
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"^[^^^ (xi < X < X2) = ^ (E {X2, x) (gi, k2, q2, h) {q2 + ih) (1 + cos 9) 

+^ (gi, k2, q2, h) {q2 + iki) (1 - 0086*)] 
+E^^ (x2,x) (gi, A;2,g2, ^4) (^2 - ^/^a) (1 + cos6') 
(gi, k2,q2, k'i) (ga - 2/^4) (1 - cos6')]) 

*P (xi < X < X2) = P^\hr\ih - A;4)sin^ [^(xi,x) (gi - 2^:2) + ^"^ (xi,x) (gi + iA:2 



den y g(x) 
den 



*P(X2 <x< 00) = --^v^M^[e^''^"""=^^(gi,fc2,g2,fc4)(l + cos^) 

cten '- 

+e'>^^(--^)^ (q,,k2,q2,ks) {I - cos9)] 
^ (X2 < X < cx)) = -^Vhq^2 [e^'^^^^-^^)^ (gi, ^2, g2, A;4) - e^'^^^^"^^)^ (gi, A;2, g2, A^a)] sin^ 

Electrons initially in the right electrode have the following wavefunctions along the 
structure : 



R 



— 00 < X < Xi , 



^T(T) 

R ' - ^ den 



:^v^g^g^^(gi,/c2,g2,/c4)cos^e ^'^^(^ 

den 2 



i{T), 



-00 < X < Xi) 



8z 



v/gig2A;3*(gi, ^i, ^2, ^4) sin _e"*^-2(^-^i) 



'^^r\xi < X < X2) 

< a: < X2) 



'4~\ ~T^'^(9i, h, g2, A;4) cos^ [-E(xi, x)(gi - iki) + E ^(xi,x)(gi + ik^ 
den V q[x) 2 



■^W^^^(gi,A;i,g2,fc4)sin^ [-E(xi,x)(gi - ik2) + E ^(xi,x)(gi + z/c2)] 



\l/]j^^^(x2 < X < 00) = cos ^ 



e 1 



\l>;[j^^(x2 < X < 00) = sin ^ 



2v^ 
6 1 



-iki(x-X2) „iki{x-X2) 



, g sin e 8gig2fc3(fci - fca) ^.fc.^-^,) 



2x/E 



den 



2^/fe 



-2A;3(x-X2) ] ik-i{;x-X2) 



cos -- 



61 sin6' 8gig2/c3(^i - ^2) 



^ifc4(a;— 2:2) 



2x/E 



den 



To obtain v|/-L(i) and "^^^^"^ from v[/T{T) g^j^^j v|>-l(T)^ q must be replaced by —9 and /ci 
(A;3) by k2 (^4) in the above formulae. 
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